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Abstract 

There is described a spacetime formulation of both nonrelativistic and 
relativistic elasticity. Specific attention is devoted to the causal structure 
of the theories and the availability of local existence theorems for the 
initial-value problem. Much of the presented material is based on joint 
work of B.G.Schmidt and the author (in Class. Quantum Grav.20 (2003), 
889-904). 

1 Introduction 

The ancient field theory of continuum mechanics, created by the math- 
ematicians J.Bernoulli, Euler and Cauchy, has grown into a subject of 
great importance both for its mathematical interest and its applications 
in material science. In its relativistic guise this theory has not been de- 
veloped very far yet. Important references are [3], |17| and 
|10| . The book |15| is also an excellent source. 

We start, in Section |21 by describing the nonrelativistic theory in a 
framework akin to that used in relativity, namely that of Galilean space- 
times. In the following Sectionj^we describe the concept of hyperbolicity 
appropriate for the resulting system of 2nd-order partial differential equa- 
tions. We then describe a way of rewriting such a system in symmetric 
first-order form. The condition of symmetric hyperbolicity for this lat- 
ter system, although true in many cases of physical interest, is however 
a more stringent requirement than that of hyperbolicity for the original 
second-order system. In Section |1] we write down a class of states, the so- 
called "natural states", which satisfy the assumptions of (both the first- 
and second-order) hyperbolic theory outlined in Section |3 so that there 
is available, for initial states sufficiently close to natural ones, a local 
existence theorem for the Cauchy problem.^ Finally, in Section |S] we de- 



^ There are in fact theorems even on global existence (see |14|') 
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scribe the necessary changes as one goes from the nonrelativistic theory 
(or rather: "GaUlean relativity") to Einsteinian relativity. 

2 Nonrelativistic Theory 

We start with a Galilean spacetime M(see e.g 0). This is furnished by 
R'*, endowed with a symmetric, degenerate contravariant metric /i^" of 
signature (0 + ++) and a choice of covector field satisfying 

/i^V, = (1) 

We also assume we are given a flat connection which annhilates both 
ft^" and Tfj,. The matter flow is described by a vector fleld v'^ normalized 
by v'^T^ — 1, i.e. of the form 

v''d^. = dt + v'd,, (2) 

where v' = v^{t,x-') and (a;'') = (t,x^) are flat coordinates in which 
h'^"d^,di, — S^^didj and Ti_idx'^ = dt. The given flat connection V is not the 
only one annihilating (h'^'^Tp). One can use this freedom to describe the 
effect of gravity by using the connection V, where V^u^ — V^lOi, — C^^^ujx 
with C^^ = T^T,yh''^'^V aU , U being the gravitational potential. One easily 
checks that the connection V again annihilates (/i'"',t^). Having de- 
scribed the kinematical arena, we now turn to the speciflc class of physi- 
cal models we consider. We start by writing down a "stress-mass" tensor. 
This tensor is not such a natural object in the nonrelativistic theory as 
the stress-energy tensor is in relativity, for two reasons: the lack of a 
non-degenerate spacetime metric and the fact that the Lagrangian of the 
nonrelativistic theory breaks the Galilean invariance. We will nonetheless 
use the concept of stress-mass here, since it greatly facilitates the task of 
moving back and forth between the Galilean and the Einsteinian theory. 
The mass-stress tensor has the form 

T"-" ^ pv^v-" +t'"', (3) 

where the Cauchy stress tensor t'^'^ is purely spatial in the sense that 
t'^'^Ti, = 0. Furthermore we have that p — nmo, where n > is the particle 
number density and mo the mass per particle. The continuity equation 
in our language takes the following form. Take e, the volume form on 
M defined in the adapted coordinates as Sijko = ^ijk, and consider the 
three- form A'^ given by A^^j^a ~ ne^.vXpv'' . Then conservation of mass is 
given simply by dN — 0. In the standard coordinates, A'^ is given by 

= neijk{dx' - v'dt){dx^ - v^dt){dx'' - v'^dt) (4) 

and dN = is of course equivalent to 

d^ijiv'^) = dtn + di{nv'') = 0. (5) 

The matter fleld equations are given by 



(6) 
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in the absence of gravity, otherwise the derivative V^j has to be used in 
Eq.lO- We immediately see that the equation r^V^T^" = is already 
implied by the continuity law Eq.Q. In order to turn Eq.Q into proper 
field equations we have to specify the dependent variables. These are fur- 
nished by maps / sending points of spacetime M into a manifold B called 
body or material manifold. This material manifold should be viewed as 
an abstract set of labels which parametrize the particles making up the 
continuum. Thus / is the "back-to-labels-map" . If we choose coordi- 
nates (X-*) with A = 1,2, 3 on B, we can write = f^{t,x^). The 
relationship between the map and the vector field is given by 

v^d^f" = dtf" + vd.f^ = (7) 

Suppose that B is endowed with a volume form ^Iabc and set 

fn = N, (8) 

where /* denotes pull back under the map /. Eq.Q defines n in terms of 
f^, namely there holds n = det{dif^). We assume the map to be such 
that f{t, .) is a diffeomorphism onto its image in B for all t and oriented so 
that n is positive. This implies that the map / is of maximal rank. Hence, 
given /, Eq.l|7| has a unique solution v^. Here, and in what follows, the 
spacetime field «^ is always viewed as a function of {f^,d^f^). We now 
introduce the concept of "strain" by means of quantities H^^ defined by 

H^^ = h^-{d^f^)id.f^) (9) 

Clearly H^^ is positive definite. Consequently there exists the inverse 
Hab defined by 

H^'^Hbc = <5^c (10) 
We now assume that the Cauchy stress tensor in Eq.JSJ is of the form 

= nTAB{dpf^){dafWh'''' (11) 

with 

-^^=2^^ (12) 

for some function e = e{f^(x), H^'~^ (x)), called stored-energy function in 
the elastic literature. As an example take the case where e just depends 
on n. Note this makes sense since n can be written in terms of H^^, 
namely 

6n^ = H^"^' H'"'' H"^"^' Q.ABC^A'B'C' (13) 
When e depends only on n one finds that 

t'"' ^ ph"-' , (14) 

where p is defined by 

p = n- (15) 

The field equations here are the Euler equations for a perfect fluid. 

For completeness we outline the proof of Eg. 1141 . We first claim that 

-^Hab (16) 
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The trick how to obtain Eg. 11161 is to first compute qJ"b H^^ , using 
Eq.l^Hl- One quickly finds that 



We next define a quantity a by 

= {d.f)h'"'HAB (18) 

One checks that 

{d^.f^)F^B = 5^B (19) 
Using Ff^AT^ = and {d^f^)v'' = 0, it follows that 

F%(a./^)=<5^-^'V. (20) 

Using Eq.'s lltil^m in Ea. lllH . we immediately obtain Ea. (ll4l l together 
with HSU . 

Let us return to the case of a general elastic solid. The field equations 
are of the following form 

tttt: — 777 dxdi. f — lower-order derivatives of / , (21) 

which, by the remark following Eq.JSJ, are equivalent to 

M^-g d^d^f = Ga (22) 

where M^^, defined as 

= i9^f)HcA^g-j^, (23) 

and Qa are functions of {f^, 9^/^). Remarkably the quantities Af^g turn 
out to satisfy the symmetry 

= M^^ (24) 

The symmetry Ea. l241 is no accident. It is due to the fact that the 
field equations are the Euler-Lagrange equations of a variational principle. 
Explicitly we find 

MTb = -Pv^v'Hab + n[TABHcD + 2tc^aHb)d + ^^^WfdJ^, 

(25) 

where we have defined d'^f^ — h'^'^d^f^. The associated Lagrangian, 
which in particular satisfies 

is given by 

C = n(—mQV^v-' Sij — e). 
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Note that the quantities n, e, should be all regarded as functions of 

This ends our description of nonrelativistic elasticity in its spatial (or 
rather:" spacetime") form. We remark that all standard treatments in the 
literature (see e.g. or 113) 1 prefer the material form based on the map 
F\t,X'^) defined by 

f'^{t,F\t,X'')) = X"^. (27) 

From the relativistic point of view the spacetime form is preferable, 
since, in a relativistic spacetime, there does not exist the standard t = 
const-foliation available in Galilean spacetime. 

3 Some Hyperbolic Theory 

An equation of the form of 12211 will be called hyperbolic if Af^^ satisfies 
the symmetry l{'24^ and the following holds: firstly there should exist a 
subcharacteristic covector, i.e. a covector so that 

-'^^abCmCi' is negative definite (28) 

Secondly there should exist a timelike vector, i.e. a vector so that 

M'^'^rjf^m^rj^m^ is positive definite (29) 

for all nonzero 77^ with X^r^f^ = 0. Let us pause for a moment to explain 
by means of an example simpler than elasticity how the change of sign 
between Ea. l281 and Ea. l291 arises. The example is the equation for wave 
maps, where A/^^ = q^^Gab with t/'"^ the spacetime metric and Gab 
the Riemannian metric on the target space. Now the notion of timelike 
has its standard Lorentzian meaning, and subcharacteristic covectors are 
timelike covectors. The sign change between 1281 and 12911 is then simply 
due to the fact that (co-)vectors orthogonal to a timelike (co-)vector are 
spacelike. 

A vector X'' is called causal if X**^^ 7^ O for all subcharacteristic 
covectors Clearly all timelike vectors are causal. In general there will 
be a gap between timelike and noncausal vectors due for example to the 
existence of different characteristic (sound) cones. A covector A;^ is called 
characteristic if the symbol of the PD operator on the l.h. side of Eg. 12211 . 
namely the quadratic form MAB{k) given by 

AfAs(fc) = Af^^fc^fc. (30) 

is degenerate. A vector X'^ is called (bi-)characteristic if it is of the form 

X" = A^^^fc.m^m^ (31) 

for a characteristic covector fc,i and such that M'^gk^k^m^ — 0. 
This characteristic vector is tangent to the characteristic sheet to which 
k belongs where this sheet is a regular surface. 
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The above definitions are essentially taken from the book [^-^ The 
notion of characteristic vector from as above is new, the classical one 
breaking down when the set of characteristic covectors (also called "nor- 
mal cone" or "slowness cone" in the literature) has singularities due to in- 
tersections between different sheets (there are in general three such sheets 
which are given by the zero- level set of the eigenvalues of Mas). Even 
in the absence of singularities of the normal cone the set of characteristic 
vectors (also called "ray cone" or "wave cone" in the literature) will in 
general have cusps (see e.g. Chapter VI of W). 

A yet more general notion of hyperbolicity, due to Kreiss, of which the 
above is a special case, is that of strong hyperbolicity (see ED)- 

A key point regarding these definitions is the availability of an ex- 
istence theorem independently of the form of the lower-order terms in 
Eq.^.^ Namely, suppose there is a hypersurface S of spacetime, to- 
gether with initial data for and on S so that, for these data, 
the surface S has everywhere subcharacteristic conormal. Then choose 
a vector field which is timelike on S, whence transversal to S. Use 
this vector field to Lie-drag S into the future. Since the properties of be- 
ing subcharacteristic and timelike are "open" conditions, we thus obtain 
a spacetime neighbourhood N G M of S, which is foliated by surfaces 
which are subcharacteristic for all maps / close to the initial one and 
where X'^ is timelike with respect to such configurations. Now take coor- 
dinates {y^jU') so that the leaves of the foliations are given by y'^ = const 
and the vector field is given by X'^d^ ~ do. In these coordinates the 
equation (I22I I takes the form 

[M'i°B^o + {M'i'B+^^^B)^o^^+^^/g^^^,]f'^ = lower-order derivatives of/-^. 

(32) 

Furthermore there should hold 

MABm^m'^ < 0, Ar/ghljm'^m'^ > (33) 

with li, both nonzero and all values of (/^, 9^/^) close to those cor- 
responding to the initial data. If the neighbourhood A'^ is of the form 
{j/" G [0,r]} X {y € R^} and the initial data satisfy some decay prop- 
erties for large \y\ one can now appeal to a basic theorem in [3| to infer 
existence of a unique solution for sufficiently small T. (We do not spell 
out the precise differentiability requirements.) The asymptotic conditions 
imposed in the above theorem are of course not always appropriate, and 
one would in any case like a local statement amounting to uniqueness in 
the "domain of dependence" of initial data in open subsets of S. Such a 
theorem is proved in 5 , the appropriate notion of domain of dependence 
being based on causal curves (in the sense of causal vectors as described 
above). Consequently the nonlocal nature of the uniqueness part of the 
theorem in |5| is in fact irrelevant. 

In PP the autors chose to cast the equations of elasticity theory into 
that of a first-order symmetric hyperbolic system, which goes as follows: 

■^We have only added the word "timelike" for vectors having the property in Ea. l29i and 
"subcharacteristic" for covectors satisfying Ea. l28i . 
^In |1] a result to that extent is stated without proof. 
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Define 5-index quantities W'^ab'-^^ by 

W'^AB^''^ := X'^M^-'ab - 2X"'M''l''flA, (34) 
where X is a timelike vector. We now replace 12211 by the following first- 
order system: 

W^-'ab^^HJ, F)dxF''. = X^if, F)gAif, F) (35) 

-X\f,F)d^f'' = -X\f,F)F^x (36) 
together with the constraint 9,,/'^ — F^^. Since 

W'^ab'^^^ =W''''b/^\ (37) 

the system I|85I36|I is symmetric. One then finds that H35I86|I is equivalent 
to the original second-order system 1221 Next recall that the symmetric 
system is called symmetric hyperbolic if there exists a subcharacteristic 
covector i.e. one so that the quadratic form defined by the left-hand 
side of the system Eg. 13511 . i.e. W^'^ ab^^^ ^\ is negative definite in the 
variables m^^. Suppose ^ is subcharacteristic for the second-order system 
and X'^^^ > 0: is it then subcharacteristic also for the system I|85I36^ ? 
The answer in general is "no" as we will see in the next section. 

4 Natural States 

We now further specify the "equation of state" given by the stored-energy 
function, as follows: We assume B to be endowed with a fiat Riemannian 
metric Gab and that the volume form Q is compatible with Gab- The 
stored-energy function e is then assumed to be of the form e = e{H^^), 



where H^^ now refers to coordinates X^ on B in which Gab ~ Sab- 
(Note this implies that the field equations 1221 and 1351 have Qa = 0.) 
More specifically we suppose e to satisfy 

e = \eabcd{H'''' - 5'*^)(jy'^° - 5'='°) + 0{{H - Sf) (38) 



for certain constants Eabcd = E(^ab){cd) = Ecdab- Clearly there are 
21 independent such constants available. One furthermore assumes these 
constants to be such that 

EABCol^m^l^m^ > (39) 

O 

for m^, both nonzero. The definition, then, of a natural state / is that 

o 

of a map which corresponds to a configuration of zero strain in that 

H^^ix) = {dj^ix)){d,/fix))h^"{x) = 5-"^. (40) 
Note this relation implies 

Coy = (41) 

*In (0 we used a special property of elasticity to prove this equivalence. It is not hard to 
see that this equivalence works generally for the system 1221 when is timelike. 
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i.e. that u is a rigid motion. The absence of linear terms in Eg. 13811 
furthermore implies that natural states are stressfree, i.e. have tab = 0. 

o 

Using also that n= 1 we find that 

M'X'b = -mo v-'Sab + EacbdR^'^R'"' (42) 

with R'^'^ = R^bS^'^ and R'^b a (in general time dependent) rotation ma- 
trix. The covector r^, from Ea. l42L is clearly subcharacteristic. Further- 
more the vector n ^ is timelike iff the conditions 13911 are valid. Therefore 
the second-order equations are hyperbolic at natural states. 

We now turn to hyperbolicity of the first-order system. Taking 
X'' =5'' and using Eq.s 142I34L we see that 

W'^-'ab^^^tx = -mo > v-Sab - EacbdR'''' R''" (43) 

Thus Tfj, is subcharacteristic for natural states iff 

EABCDm'^''m'^'' > (44) 

for all nonzero elements m^^ = m^^^\ which is clearly a stronger re- 
quirement than 1391 . 

The only case which is easy to analyze fully is the isotropic case ^ 
where 

Eabcd = >^5ab5cd + 2^5c{a5b)d (45) 

The constants A,/i in Eg. 1451 are the standard Lame constants. They 
should not be confused with indices (p,v). The "rank-one convexity" 
condition Eg. 13911 is equivalent to 

moC2 = /i > 0, moCi = 2/i + A > (46) 

The eigenvalues of -^MAB{k) relative to Sab, which are real of course 
since Mab is symmetric, are given by 

Ai(fc) = cf ff^^fc^fc, (47) 

and 

\2{k) = A3(fc) = cl g'^-'k^K, (48) 

where 

ff"" = /i^" - v^" v\ l^"" = ftf" - 1 > °v\ (49) 

(A Lorentzian metric of the above form is nowadays called "acoustic met- 
ric" or "Unruh metric".) Thus the normal cone consists of two sheets. 
The one corresponding to Ai is associated with a longitudinal ("pres- 
sure") mode propagating at speed ci, the second one corresponding to 
two transversal ("shear") modes of speed C2. If C2 < ci, the first sheet 
lies inside the second sheet. Subcharacteristic covectors, for which both 
Ai and A2 are negative, lie inside the inner cone. One such covector is 
r^. In fact, lies on the central ray inside the two cones in the follow- 
ing sense: the vector ii'' defines a family of parallel hyperplanes in the 

^In we mistakenly interchanged the constants A and in equation (4.16) of that paper. 
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cotangent space. These intersect the normal cones in 2-surfaces which, in 
the metric h'^" , are standard spheres centered at the point where the ray 
of intersects this hyperplane. The ray cones dual to the above normal 
ones are given by the equations 

(50) 

where /i^i, is the unique tensor satisfying 

h^^h^p = p — v'^Tp, hp^^v" = (51) 

12 12 

Note that g^^, (resp. 3^^) are the inverses of <? (rasp, '"^). Clearly, 
the cone of shear waves is now the one lying inside. Timelike vectors X, 
for which all covectors k with X'^kp, = have Ai and A2 both positive, lie 
inside this inner ray cone. One such timelike vector is v'^ , in fact, it lies 
on the central ray of the two ray cones in a fashion exactly dual to that 
explained for Tp. Causal vectors may be "faster" in that they lie inside or 
on the outer ray cone. 

We now turn, finally in this section, to the question of hyperbolicity 
of the first-order theory at natural isotropic states. The condition II44I I is 
valid if and only if 

c? > ^ci > 0, (52) 

which is physically entirely reasonable, since elastic materials typically 
have Ci/c2 approximately 1,7. But one can do better than that. One first 
notices that both the symmetries of and the equation I22II remain 

untouched if the quantities M'^^ are replaced by given by 

M^/^^M%l+K^:^, (53) 

where A^"^ = ^as ~ ^[ab]-^"^ ^^^^ ^'^^^ replacement can be viewed as 
coming from adding a total divergence to the underlying Lagrangian (see 
[5|)). While it is easily seen that second-order hyperbolicity is unafi^ected 
by this replacement, this is not the case for the associated first-order 
system. In the case at hand we can, by adding to Eacbd a term of the 
form 

(4c^ - 2e)5AycSD]B, (54) 

arrange for Ea. l44ll to hold if we take e in the range < e < min{2c2, -5^), 
which of course can alway be satisfied when ci and C2 are both non-zero. 

Let us point out that the standard equations of linearized elasticity at 
a natural state are obtained by simply freezing the coefficients in Eg. 12211 
and setting the right-hand side equal to zero. The Cauchy problem for 
these equations is studied e.g. in 0. In the isotropic case one finds that 
the solution at {t,x^) does not depend on data at t = inside \x\ < C2t, 
i.e. inside the past inner ("shear") ray cone. 

5 Relativist ic Theory 

Having available the spacetime form of the nonrelativistic theory it is easy 
to write down its relativistic version, so we will be brief, merely pointing 
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out the necessary changes. We start out with a relativistic spacetime 
(M, g^i/) with g^i/ a Lorentz metric. The configurations are again maps 
from spacetime to B, the latter endowed with a Riemannian metric Gab 
and compatible volume form Habc- The maps / should be of maximal 
rank and such that the inverse image under / of each point of B in the 
image of / is a timelike curve in M. Thus there are timelike vectors u'^, 
unique up to scale, so that 

{d,f^)u^ = (55) 

We denote henceforth by u'^ the unique solution vector of Eg. 1)55^ which 
is future-pointing and normalized by g^^u'^u'^ = —1. The quantities H^^ 
are defined as 

H^^ = (56) 
and are again positive definite. The particle number density n is defined 

by 

= nA'B'c'^ABcH'^'^' H''''' H'^'^' , n > (57) 

The Lagrangian of the theory (we set the speed of light equal to one) is 
taken to be 

£ — —p — —n{mo + e), (58) 

with e a function of {f^,H^^'). Varying C with respect to g^" gives the 
stress energy tensor 

T""" = pu^u" + 2nTAB{d''f'^){d'' f), (59) 

where tab = as before and d^f^ = g'"'dr,f^. The field equations 

are again of the form Eg. 1221 with 

M^-ab = -liABU^n- + UacbdF^^F'"', (60) 

where 

PAB = pHaB + riTAB (61) 

and 

Uacbd = n(TABHcD + tacHbd + tbdHac + '^ qj^ac ) '2pHaicHd]b- 

(62) 

The last term in Eg. (1621 does not contribute to the equations of motion. 
The second term on the right in Eg. 1611 is clearly a relativistic contribu- 
tion. 

There are slight complications in curved spacetime regarding the no- 
tion of a natural state: in order for a natural state to exist, the spacetime 
metric would have to allow Born rigid motions and the material metric 
Gab would have to be isometric to the metric on the quotient of A/ by the 
action of this motion. We avoid this difficulty here by confining ourselves 
to special relativity and taking the natural configuration to be at rest in 
some inertial system. Thus we assume (M, g^i/) to be Minkowski space. 
We also suppose the "natural" motion to be of the form u'^dfj, = dt in co- 
ordinates {t,x^) in which g^^dx'^dx^ — —dt^ + Sijdx'dx-' . We also assume 
Gab ~ Sab as before. Using 1401 . a natural map corresponds to a time 
independent rotation in an inertial system which, in the isotropic case. 
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can be taken to be the identity without loss. One then merely replaces, 
in the expressions I|47I48^ for the different cones, the covector r,j by the 
covector ~ u^— —g^i, u'^ and the symmetric tensor h'^" by h^" given by 

/jM-- ^^M-_^ °M (63) 

Furthermore one replaces the tensor /i^i^ in H51|l by hfj,v given by 

h^v = g^v+ Ufj,u„ . (64) 

Thus, writing c for either ci or C2, the associated normal cone is now given 
by the Lorentz metric 

^M.^gM. + (l_^) °M °. (65) 

and the ray cone by the inverse metric, namely 

g^iu = fffii' + (1 ^ c^) "m"!' • (66) 

One easily infers from these relations that in the present coordinates the 
special relativistic equations, linearized at a natural state, are exactly 
identical with the nonrelativistic ones, when the latter are written in co- 
ordinates where = dt with v'^ an inertial motion. Of course the 
geometrical objects in both theories are different - and thus the behaviour 
of the two linearized theories under change of coordinates is also different. 
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